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Abstract 

The use of atomic decay rates as a probe of sub-vacuum phenomena will be studied. Because 
electromagnetic vacuum fluctuations are essential for radiative decay of excited atomic states, decay 
rates can serve as a measure of the suppression of vacuum fluctuation in non-classical states, such 
as squeezed vacuum states. In such states the renormalized expectation value of the square of 
the electric field or the energy density can be periodically negative, representing suppression of 
vacuum fluctuations. We explore the extent to which atomic decays can be used to measure the 
mean squared electric field or energy density. We consider a scheme in which atoms in an excited 
state transit a closed cavity whose lowest mode contains photons in a non-classical state. The 
change in the decay probability of the atom in the cavity due to the non-classical state can, under 
certain circumstances, serve as a measure of the mean squared electric field or energy density in the 
cavity. We derive a quantum inequality bound on the decrease in this probability. We also show 
that the decrease in decay rate can sometimes be a measure of negative energy density or negative 
squared electric field. We make some estimates of the magnitude of this effect, which indicate that 
an experimental test might be possible. 
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I. INTRODUCTION 



It has been known for some time that negative energy densities and fluxes are a generic 
prediction of quantum field theory States which involve negative energy, such as the 
Casimir vacuum state and squeezed states of light, have even been produced in the labo- 
ratory 0, [3[ . Negative energy density may be viewed as an example of a sub- vacuum phe- 
nomenon, whereby the vacuum fluctuations are suppressed below their level in the Minkowski 
vacuum state. The experiments which have been done measure an indirect effect, such as a 
Casimir force or a change in photon counting statistics, not the energy density itself. This 
raises the question of whether more direct detection of negative energy density or related 
effects is possible. 

The gravitational effects are far too small to be feasible in a laboratory experiment. 
Furthermore, the magnitude and duration of negative energy densities and fluxes are con- 
strained by quantum inequalities @, @, @, 0, II] • These are constraints derivable directly from 
quantum field theory, which yield an inverse relation between the magnitude of the negative 
energy density and its duration. Marecki [9| proved quantum inequality-type bounds on 
the magnitude and duration of squeezing in quantum optics experiments involving squeezed 
states of light. These constraints make the search for local measurements of sub-vacuum 
effects more challenging. 

Nonetheless, there seems to be no barrier in principle to constructing a negative energy 
detector which relies solely upon electromagnetic interactions. One model for such a detector 
was proposed by Ford, Grove, and Ottewill [10], who analyzed a system of atomic spins, 
placed in an external magnetic field and coupled to a quantized electromagnetic field. When 
the system is coupled to a non-classical state of photons, such as a squeezed vacuum state, 
there can be a transient increase in the average magnetic moment, as compared to when the 
quantized field is in the vacuum state. This can be viewed as a suppression of a depolarizing 
effect of vacuum fluctuations on the spins, resulting in momentary "re-polarization". In 
general, this change is not directly correlated with energy density. However, for a certain 
ranges of parameters of the system, the change in magnetic moment is in phase with the 
periods of negative energy density. Under these circumstances, this spin system represents 
a non-gravitational negative energy detector. Unfortunately, the fractional change in the 
mean magnetic moment is quite small unless the photon energies approach the 7-ray range. 
Consequently, it is questionable that one could find a way to detect the extremely rapid 
changes in the mean magnetic moment. 

Other models have been treated by Davies and Ottewill [ll|, who studied the detection of 
negative energy fluxes using various types of switched monopole particle detectors, building 
on earlier work by Grove Marecki and Szpak [l3[ modeled spontaneous light emission 
from two-level atoms coupled to a quantized electromagnetic field. They derived a Volterra- 
type equation which controls the time evolution of the amplitude of the excited state. 

In this paper, we follow the spirit of these earlier works and look for an indirect "tracker" 
which might lead to the detection of negative energy or vacuum fluctuation suppression- 
type effects. It is well-known that electromagnetic vacuum fluctuations are essential for 
the spontaneous decay of excited states of atoms. Without the coupling to the quantized 
radiation field, all energy levels of an atom would be eigenstates of the Hamiltonian, and 
hence stable. This suggests that a suppression of the usual vacuum fluctuations could 
increase the lifetime of an excited state in a way which might be observable. 

We consider a model of a two-level excited atom interacting with a quantized electro- 
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magnetic field, using first-order perturbation theory. The field is confined to a closed cavity 
with two of its dimensions much larger than the other one. This is so that we can arrange 
the transit time of the atom through the cavity to be small compared to the light travel 
time across the larger two dimensions of the cavity. Our idea is to prepare the cavity field 
in a non-classical state, such as a squeezed vacuum, and fire the excited atom through the 
cavity along its shortest dimension, before the state of the field has time to change very 
much. We want the atom to interact with the field during the period when the renormalized 
expectation value of the square of the electric field, (E 2 ), is negative. The purpose is to 
see whether this will suppress the de-excitation probability of the atom, compared to its 
vacuum value, in a way that will be correlated with the periods of (E 2 ) < 0. If so, our 
system would function as a negative E 2 detector. We work primarily with cavity modes, 
for which (E 2 ) ^ (B 2 ), so periods of (E 2 ) < do not generally correspond to periods of 
negative energy density. However, under certain conditions, our model will also serve to 
measure the energy density. 

For most treatments of this type in quantum optics, the Jaynes-Cummings model for the 
interaction between the atom and the field is used. This model makes use of the rotating wave 
approximation, which effectively ignores terms with rapidly oscillating exponentials [l4( . We 
explicitly do not make this approximation, because the effects we are interested in are highly 
transient ones, and depend on these terms. 

It also should be mentioned that any scheme which is designed so that a particle or 
observer interacts only with negative energy in flat spacetime can be ruled out by the 
quantum inequalities and the averaged weak energy condition, which is known to hold in 
Minkowski spacetime. In a realistic version of our atom-cavity system, the atom would 
have to pass through holes cut in the cavity walls. Even in the event that the atom passed 
through net negative energy while passing through the cavity, edge effects from the holes 
will contribute enough positive energy to satisfy the averaged weak energy condition. This 



effect has been discussed recently by Graham and Olum [15|, and by Fewster, Olum, and 



Pfenning [16 



In our case, we are comparing the de-excitation probability of the atom in an excited state 
of the cavity to when there is only vacuum in the cavity. In particular, we are interested in 
situations where this probability is suppressed relative to its vacuum value, and where the 
periods of suppression are in phase with periods of (E 2 ) < 0. We will be concerned with 
changes in (E 2 ) or the energy density due to changes in the quantum state in the cavity. 
In this case, edge effects due to the holes in the cavity will cancel out. It is possible for 
the difference in the net energy seen by an observer to be negative, as was shown by Borde, 
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Throughout this paper, we will regard (E 2 ) and the energy density as being set to zero in 
the vacuum state of the cavity. Thus we are concerned only with changes due to changing 
the quantum state of the cavity, and not with Casimir-type effects due to the geometry of 
the cavity. It is well known that the presence of boundaries can make the energy density 
or (E 2 ) smaller than in empty space vacuum state. However, this effect on decay rates of 
atoms is difficult to distinguish from the effects of the changes in mode structure, such as 
the change from a continuous to a discrete spectrum. The effects of cavity geometry on 
atomic decays has been extensively studied in recent decades, with an e arly treatment given 



by Babiker and Barton [18j • For a recent review, see for example Ref . [19[ . 

The outline of this paper is as follows: In Sect. [Til we develop some formalism for de- 
scribing the interaction of cavity modes with an atom. This is further developed in Sect. IHfl 
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where we obtain detailed expressions for the de-excitation probability for an atom travers- 
ing a cavity with one mode excited. In Sect. IIVI we make some numerical estimates of the 
size of the effect and discuss the feasibility of observing it. Our results are summarized in 
Sect. El In the Appendix, we derive a quantum inequality- type bound on the de-excitation 
probability. Unless stated otherwise, we work in Lorentz-Heaviside units where % = c = 1. 



II. ATOM-CAVITY INTERACTION 

Consider a two-level atom with states and \ip 2 )- If the atom passes through a cavity 
containing a quantized radiation field, then the final state of the system will in general be 
an entangled one. First, consider the case where the initial photon state in the cavity is 
a single-mode number eigenstate, = \n), and the atom is in the state \4>i). Then the 
initial state of the system is 

\%) = In) IVi) ■ (1) 
If we consider the case where the interaction changes the photon number by at most one, 
then the most general final state is of the entangled form: 

\V f ) = B x |n) |Vi)+S 2 |n+l) \ii 2 )+B 3 \n-l) \ipi}+B 4 \n+l) |Vi>+£ S \n) \ip 2 )+B 6 \n-l) |^ 2 ) . 

(2) 

Suppose we want the probability of finding the atom in state \ip 2 ), irrespective of the photon 
state. Then we should project |\&/) onto the subspace where \ip a tom) — \^2)'- 

(H*f) = B 2 \n + l) + B 5 |n) + B 6 \n - 1) . (3) 

However, B 5 = 0, because the interaction term only connects photon number states which 
differ by ±1. 

The interaction Hamiltonian for the atom and the field in the cavity, in the dipole ap- 
proximation, is given by: 

H' = — d • E s (4) 

where d is the dipole moment of the atom and Eg is the Schrodinger picture electric field 
operator for the quantized cavity field, evaluated at the atom's position. The dipole ap- 
proximation assumes that we can ignore the variation of the field across the size of the 
atom. The Schrodinger picture electric field is related to the Heisenberg field E(x, t) by 
Es(x) = E(x, 0). The Heisenberg field has the mode expansion 

E(x, t) = £K A e fcA / feA (x) e-*"* + ahxe k x / feA (x) e*"*] . (5) 

k\ 

Here fk\{^) is the spatial part of the mode function, which we take to be real, and e^x is a 
linear polarization vector. 

First order perturbation theory yields 

B 2 = -i r dt' (n+ l,i) 2 \H'\n,ipi) e 1 ^-^' 

Jt 

= iV^+l f 1 dt'e kX ■ (H^i) /*a(x(0) e l ^- A£)t ' (6) 

Jto 

B 6 = -i t dt'in- l,ip 2 \H'\n,i>i)e- i{ui+A£)t ' 

= i^i r dt'e kX ■ (^|d|V»i> /*x(x(0) e-'^ +A ^' . (7) 

Jto 
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We assume that \ipi) is the higher of the two energy states, and let Ae > be the energy 
difference between the two atomic states. The probability of finding the atom in ^2) is 

P 2 = \(H^f)\ 2 = \B2\ 2 + \B 6 \ 2 . (8) 

Note that the contribution of B G would be ignored in the rotating wave approximation. 
Let us now consider the case of a general one-mode initial photon state 

00 

H) = E C « |n> • (9) 

n=0 

The most general final state of the atom-cavity system is 

l*/> = E l^l) + A 2n\n) m) ■ (10) 

n=0 

Now the projection of \^ff) onto l^) is the vector in the photon state space given by 

00 

y> 2 \* f ) = y EA 2n \n), (11) 

n=0 

and the probability of finding \ip 2 ) is 

00 

p a = l<v*l*/>l 2 = E \ A ^\ 2 - ( 12 ) 

n=0 

Let T be the transition matrix, so that Tj, = (f\T\i) is the amplitude to make a transition 
from state \i) to state |/). Then we have that 

00 

= ^ c n (mi})2\T\nil)i) . (13) 

n=0 

This latter matrix element between energy eigenstates is then given in first-order perturba- 
tion theory by an integral of a matrix element of H' . 



III. DE-EXCITATION PROBABILITY FOR AN ATOM IN A CAVITY 

Consider an atom passing through a rectangular cavity with dimensions a, b, and d aligned 
along the x, y, z axes respectively, as shown in Fig. [TJ We will assume that b < a < d, with 
b ^ d, and that the velocity of the atom is parallel to the 6-dimension of the cavity, v = vy. 
We will also assume that v <C 1, so that we may ignore relativistic effects. 



A. Perturbation Theory Results 

Using first-order perturbation theory, we can write the transition matrix element between 
two states of definite photon number as 

(fj#2|T|7#i) = -i [ h dt' (mip 2 \H'\n'ip 1 ) e iAEayst ' , (14) 

Jt 
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FIG. 1: A cavity for standing wave modes. We require that b < a < d, and b <C d. The atom is 
assumed to pass though the cavity in the y-direction, as illustrated. 

where AE sys is the change in the energy of the atom-cavity system. The matrix element for 
the interaction Hamiltonian is 



(mip 2 \H'\mpi) = - (ip 2 \dy\ipi) (m\E y \n) 

= -(ip2\dy\ih) (fVn5 m ,n-i + fVn + 1 S m ,n+l) ■ 

Here we assume that the electric field in the cavity is polarized in the ^/-direction. 
The transition matrix element, Eq. (|14p can then be written as 



(15) 



{mip 2 \T\n^i) = i (V^KI^iH/v^m.n-i^i + fVn + 1 5 m , n+ i J 2 



where 



and 



h 



rt-i P — i(u>+As)ti _ „—i(u}+Ae)to 



'to 



fh dt'e^- A ^' = - 



-i{uj + Ae) 



'to 



(16) 
(17) 
(18) 



i(u - Ae) 

Our Eq. (Tl3l) . assuming a general one- mode initial photon state, Eq. (jHJ), becomes 

-42m = i(^KlV'i) [Vm + 1 c m+ i //i + y/m c m _i /J 2 ] . (19) 
The probability, P 2 of finding the atom in the lower energy state l^) is 

P2 = Y.\ A 



l2m 



K^KlV'i)! 2 \cm+i\ 2 \f\ 2 \h\ 2 + m \cm-x\ 2 f\h\ 2 

m=0 

+2Re(y/m(m + 1) f 

I\ 1% c *m+l c m-l 



(20) 
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We use the relations 



5^(m+ 1) |c m+ i| 



(n) 



^ m |c m _i| 2 = (n) + 1 



E 



mm 



+ 1) c m _i c" 

m+l 

= E t/(n+l)(n + 2) C n C n +2 



(21) 



Our expression for the probability is then 

P 2 = \(Hdy\M 2 [(n)f 2 \h\ 2 + ((n) + l)f 2 \I 2 



+2Re(^r ^/( n + l){n + 2) c n c\ +2 f 2 h*h 



n=0 



(22) 



In the case when the initial state of the field is the vacuum state, [ji) = |0), we have 

P 2 = P 2 (U) = \{^ 2 \d y \^ x )\ 2 f\h? . (23) 

It is important to note that this form of ^(0) assumes that the atom decays into only one 
mode. This should be a good approximation near resonance, but not otherwise. Let us now 
consider the ratio 



Pi 



(n> 



+ (n) + 1 + 2 \/(n + l)(n + 2) Re 



n=0 



(24) 



From Eqs. (TTTj) and (U8i) . we have that 



1 - cos[(cj + Ae)(ti - t )] {u-Aef 



1 - cos[(cj - As) (ti - t )} (uj + As) 



2 ' 



(25) 



In the limit when the transit time of the atom through the cavity is very short, t\ — ► to, and 
we have 

'"> 

h 



1 



(26) 



In the limit when As —>■ uj, 



hit 



1 — cos 2uj{t\ — to) 
2u; 2 (t 1 -to) 2 

I (tx - to) ( 

~2 u ^ 

(ti-t ) 2 , 



and therefore we have that 



hit 



i e 



2to (ti -t ) 



3 2io;(ti-to) 



(27) 



(28) 
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In addition to letting Ae — > u, if we now also let t\ — > to, the last equation above reduces to 



hit 



i e 



2iojto 



-2iuj(t\ — t ) = e 



\I 2 \ 2 " 2w(ti-t ) 
Thus if we assume (ti — t )u <C 1 and (ti — t )Ae <C 1, then 



(29) 



1 . 



(30) 



If in addition, we take Ae m u, we have 



hit 

r 

\T I 2 " 
\ 1 2\ 

Therefore, the ratio of probabilities becomes 

P - = ■>{„) + ] 4 2 

n=0 



(31) 



^2(0) 



2(n> + 1 + 2 £ v /(n + l)(n + 2) Re[c n c* n+2 e 2i ^°) . 



(32) 



Let us compare this with the expression for the expectation value, (E 2 ), of the normal- 
ordered squared electric field operator in the cavity. We can calculate this using the ex- 
pansion of the operator in Eq. (jSJ) and the state given in Eq. for a single mode. We 
have 

(£ 2 (x ,t)> = (7 t |£ 2 (x ,t)|7*> = / 2 (x ) £ c* n c l (n\2a^a + a 2 e- 2 ^ t + {<Jf e^l) 

nl 



/ 2 (x ) c *nCi (2 n6i n + J(n + 2) (n+ l)^ n+2 e 



,— 2iorf 



+ Jn(n - l)5i, n _ 2 e 



2ioj< 



/ 2 (x ) £ [2^|c n | 2 + \j{n + 2)(n + 1) c* n c n+2 e" 



2iuit 



+ ^n(n - 1) c „c n _2 e 
If we relabel n — > n + 2, we can rewrite the last sum in the following way: 



(33) 



£ Jn(n - 1) c* n c n _ 2 e 



„2id;t 



n=2 



£ ^/(n + 2)(n + l) c* n+2 c n e 2 

n=0 

£ ^/(n + 2)(n + l) (c„ +2 c* n e~ 2 ^) 

n=0 



Using this and the fact that 

we may write Eq. (13"3"|) as 

(£ 2 (x ,t)) = / 2 (x ) 



£ n | Cn | 2 = ( n ) 



n=0 



2 (n) + 2 £ J(n + 2) (n + 1) Re(c n c* n+2 e 



(34) 
(35) 

(36) 



n=0 



Near resonance, i.e., in the limit As — > c<j, and in the limit of short transit times for the 
atom, i.e., t± — *> to, the ratio of the de-excitation probabilities, P 2 /P 2 (0), can therefore be 
written in terms of (P 2 (xo,t)) as 



P 



A(0) 



1 



/ 2 (x ; 



(P 2 (x ,t)>. 



(37) 



Therefore we see that ^2/^2(0) < 1 when (E 2 ) < 0. In the Appendix, we show that (E 2 ) 
is bounded from below by — / 2 (xo). This guarantees that ^2/^2(0) is non-negative, as 
required. 

Next we consider the case where we are near resonance, Ae — ► u, but with no restriction 
on ti — 1 . The ratio of probabilities is 



^2(0) 



(n) 



(I* 1 \ 
CnC\ +2 — ^ , 



(38) 



where 



h 
h 

hit 



1 — cos2o;(ti — t ) 



2u 2 {t x -t Q y 



1 e 



2co (tj - 1 ) 



3 2io;(ti-to) 



(39) 



In this case, Eq. f[38|) becomes 

P 2 / \ , 1 , / \ 1 — cos2o;(t 1 — 1 ) 



^(0) 



(n) + 1 + (n) 

00 



2o; 2 (t! - 1 ) 
^/(n + l)(n + 2) 



Pe 



ic n c* n+2 e 2 ^* ( e 2 ^ (tl -* o) - l) 



„=o w (*i - *o) 

Note that (P 2 ) contains the factor e 2luJt and the integral of this expression, 



(40) 



ti 



„2iu;t 



cZt e"~" = (e 

t 2iu 



1 - e 



2iuito "> 



(41) 



does not have a factor of l/(ti — to)- Thus in general, ^2/^2(0) does not seem to be 
proportional to either (E 2 ) or its time integral. When u(t\ — to) ^> 1, we have 



P> 



^(0) 



(n) + l + 0\ 



o;(ti - t ). 



+ {n)0\ 



,a; 2 (ti -t ) 



(42) 



Therefore, in this limit ^2/^2(0) > 1. 

The advantage of considering the ratio P 2 /P 2 (0) is that the precise form of the mode 
functions does not matter, since they cancel out. However, the drawback is that our ex- 
pression, Eq. (157)) . is only valid near resonance since we assume decay into only one mode. 
To avoid this limitation, it is also useful to look at the difference, AP 2 = P 2 — P2(0), so 
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that the contribution of the unexcited modes cancels, and hence we do not need to be near 
resonance. The difference, AP2, is given by 

AP 2 = |(^2|d w |Vi>| 2 / 2 [(n>(|i'i| 2 + l^| 2 ) 

+ 2 y/(n + l)(n + 2) Re(c n c\+ 2 I\I 2 ) , (43) 

n 

where 

2 1 - cos(^ + Ae)(ti - t ) 

Ml ^ O 1 

' ' (w + As) 2 

2 _ 1 - cos(^ - Ae)(ti - t ) 
11 (lu-Ae) 2 

e 2iu)tx _|_ e 2iut _ e i[u(ti+to)+Ae(ti—t )} _ e i[u(ti+t )— Ae(ti-t )} 

'" lh = A?-^ • < 44 > 

We have assumed that the mode functions, /, are real. However this assumption is not 
really necessary, since one can always absorb the phase of the /'s into the phases of the 
complex coefficients, i.e., the c's. 

In the case where we are far below resonance, uo <C As, 

2 

Ihf ~ \h\ 2 ~ ^2 [1 - cos Ae(t! - t )] , 

|2 _ e *Ae(ti— to) _ g-iAe(ti— to)j 

A/ 2 « ^ 

2 

~ — y[l-cosAe(*i-t )]- (45) 



This assumes that u>t and ut\ are both much less than 1. Physically, we require that 
w{ti — to) "C 1, and then for convenience set t = 0. Using the above expressions, AP 2 
becomes 

2 i2 [l - cos As{t\ - t )] 



AP 2 « 4|(^ 2 |dJVi)l / ^ — x [<n)+Ev( 7l + 1 )( 7l + 2 ) M c " c *™+ 2 

- 2|(^Wi>| 2 t 1 - 008 ^-^ (^(xo,t )). (46) 

Thus we can have AP 2 < 0, if (E 2 ) < at the time that the atom transits the cavity. In the 
Appendix we show that (E 2 ) is bounded from below and that, as a result, we have Eq. (1A9j) 

2|^ 2 K_hMI 2 / 2 (x ) 

As 2 

This result provides a limit on the degree to which sub-vacuum effects may suppress the 
decay probability. This limit is analogous to the quantum inequality bounds on negative 
energy densities and fluxes 0, H, @, 0, HJ. 

In the case where we are far above the resonant frequency, uj ^> As and when 
Ae(ti — t ) 1, we have 

N 2 ~ \h\ 2 « ^[l-cosc^-to)], 

I\h - - {C J ] , (48) 
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and 

AP 2 « ^KVaMvlV'i)! 2 !/! 2 {2(n)(l - coswfa - to)) 

- \j{n + l){n + 2) Re[c n c\ +2 (e^ - e^f)} } . (49) 

n 

Although it is possible to have AP 2 < 0, periods of AP 2 < do not seem to be correlated 
with either (E 2 ) or its time integral in this case. 

B. Specific quantum states 

In this subsection, we will discuss two specific non-classical states of the photon field. 
1. Vacuum plus two photon state 

One of the simplest examples of a quantum state which exhibits negative energy density 
and (E 2 ) < is a coherent superposition of the vacuum and a state containing two photons 
in the same mode. Such a state can be expressed as 

|7> = -Jnp(|0>+/3|2», (50) 

where we take (3 to be a real, non-negative parameter. In this state, the mean squared 
electric field may be written as 

(P 2 (x , f)> = f(x ) [2P + y/2 cos(2^)] , (51) 

which reaches its minimum value when cos(2u;i) = —1. If /3 < y/2/2, then this value is 
negative. At this point, the ratio of probabilities may be expressed as 

1 + 7~7~7w (2/? + >/2) • (52) 



P 2 (0) 1 + 

This ratio is plotted in Fig. El from which we can see that P 2 /P2(0) reaches its minimum 
value of about 0.55 when j3 ~ 0.32. Thus it is possible to achieve a 45% reduction in the 
decay probability in this state. 

The mean squared electric field is plotted as a function of time in Fig. [3] for the case 
(3 = 0.32. Note that (E 2 ) < about 1/3 of the time. However, because (E 2 ) oscillates at 
an angular frequency of 2u, the duration of the interval when (E 2 ) < is approximately 
At w n/(3u) w 1/uj. 



2. Squeezed vacuum state 

Another example of a quantum state which exhibits sub-vacuum effects is the squeezed 
vacuum state, described by the complex parameter ( = re**. This state has been studied 
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FIG. 2: The ratio of the decay probability in the vacuum plus two photon state to that in the 
vacuum state is plotted as a function of the parameter (3. 

<E 2 >/f 2 (x ) 




2 co t 



FIG. 3: The mean squared electric field in the vacuum plus two photon state is plotted as a function 
of time for the case (3 = 0.32. 



extensively, beginning with the work of Caves [20]. Using the results in Ref. [20j, one may 
show that the mean squared electric field in a squeezed vacuum state becomes 

(E 2 (xo,t)) = 2/ 2 (x ) [sinh 2 r + coshr sinhr cos(0 + 2ut)] . (53) 

This quantity is most negative when cos(0 + 2ut) = —1, at which point we have 



(E 2 (x ,t)) = -/ 2 (xo)(l-e 



-2r\ 



(54) 



If we compare this relation with Eq. (13 7p . we see that P<i — > for r ^> 1. Thus, in the limit 
of large squeeze parameter, the decay rate can momentarily go close to zero. It is also of 
interest to note that in this limit, the inequality Eq. f|4Tl) becomes an equality, as may be 



seen from Eqs. P6l) and (IA8I) . 
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C. Standing wave modes 

To make the time interval when (E 2 ) < as long as possible, we are interested in the 
lowest frequency mode of the cavity depicted in Fig. 1. A straightforward calculation using 
the formalism of Chapter 8 of Jackson [2l[ shows that, with our condition that b < a < d, 
the lowest frequency mode is the TE mode with p — I — 1, m — 0, where 



UJ = 7l \l^ + ^ ( 55 ) 



a 



and 



E x = E z — 0, 

UJa /7T \ /7T \ 

E v = — A 10 sm —x sm — z) , 

y 7r v a > ^d 1 

a /7T \ /7T \ 

B x = i- A 10 sm {-x J cos (-zj , 
B y = 0, 

£> z = — Zv4 10 cos (— ^) sin > (56) 

where the electric field is taken to be polarized in the ^/-direction. Here A w is a real nor- 
malization constant, which we will now determine. 

In Lorentz-Heaviside units, we can write the energy density for a classical electromagnetic 
field as 

p = ±(E 2 + B 2 ) = ±(E y 2 + \B x \ 2 + \B z \ 2 ) 



1 . [UJ 2 a 2 . /7T \ . /7T \ O 2 . /7T \ o/7T 



A, 2 



2 - no 



sin 2 (— sin 2 ^^) + ^- sin 2 (— cos 2 ( 



7r 2 v a 



+cos 2 (-x) sin 2 (^z) , (57) 



where we have used the expressions given in Eq. (156]) above. Let us normalize the vacuum 
mode functions by setting 

/t>a rb rd I 

d 3 x p = dx dy dz p = —u . (58) 
Jo Jo Jo 2 

This leads to the result 

abd(l + 



and 

2 2a; 3 a 2 
itiy — - — sm 



7T 



D. Difference in probabilities for low u 



We may use the above results to give expressions for the quantity AP2 in the case u <C Ae. 
Combining Eq. (|46p . and Eqs. fl60|) and fl55|) for the p = I = 1, m — 0, TE mode, we have 



AP 2 « Svr fl + gV /2 Kj^!^ ^(EJ sin2 (- 

I rf2 n 2hrl ( \ c \ 2 \n ) \d 
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x [1 - cos Aefa - t )] x [(n) + ^/(n + 1) (n + 2) i?e(c„ c* n+2 



(61) 



E. Does (£ 2 ) < imply (p) < 0? 

Ideally one would like our system to be such that (E 2 ) < when (p) < 0, so that (E 2 ) 
tracks negative energy density. For a classical single plane wave mode, E 2 = B 2 , in Lorentz- 
Heaviside units. Therefore, for the quantized field, (E 2 ) = (B 2 ), and p = ((E 2 ) + (B 2 ))/2 
imply that (E 2 ) < when (p) < 0. Unfortunately, for modes in a cavity the situation is 
more complicated, and in general (E 2 ) ^ (B 2 ). We may rewrite Eq. (jSJ) as 



E = E(^ E J 



(62) 



and the corresponding expression for the magnetic field operator as 

B = £faBj + a/Bj*). 



(63) 



Here 
and 



%/ 3 -(x) e 



-iuit 



Bj 



ejXV/,-(x) e 



(64) 



(65) 



The latter expression follows from the Maxwell equation VxE = — B, and the vector 
identity Vx(e/) = exV/ for a constant vector e. 

For the case where only a single mode j is excited, the normal-ordered expectation values 
of the squared fields are 



(E 2 ) = 2(a f a) |E/ + 2Re((a 2 ) E 2 



and 



(fi 2 ) =2(a + a) |B J f + 2 J Re((a 2 )B 



(66) 



(67) 



For a general mode function, the second terms in the expressions for (E 2 ) and (S 2 ) are not 
in phase. 

For example, consider our situation of interest, a single TE mode in a rectangular cavity, 
with p = I = 1, m = 0. For this mode we have, from Eq. (1561) . that 



E 



B 



E„ 



u 2 a 2 



7T Z 



Ai 2 sin 2 ^— x^j sin 2 (-7^) e 



2iu>t 



Bl + B = -A 
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cos I —a; sin —2; + — sm —a; cos 



a 
J 2 



'7T 

-a 



-2iu)t 



(68) 



These two expressions have opposite signs, so when (E 2 ) < 0, (B 2 ) > 0, and vice versa. 
Therefore, periods when (E 2 ) < do not, in general, necessarily correspond to periods of 
negative energy density. 

However, there are special cases when |(-E 2 )| ^> \(B 2 )\ and p ~ (E 2 ). For the mode 
discussed above, this occurs when x ~ a/2 and either z ~ d/2 or a d. In these cases, 
experiments which measure (E 2 ) are also detecting the mean energy density, p. 
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F. Comparison with the Model of Ford, Grove, and Ottewill 

In this subsection, we will compare some aspects of our model with the spin model of 
Ford, Grove, and Ottewill [ltj, which was summarized in Sect [H Ford et al assumed plane 
wave modes, so p = (E 2 ), whereas we use cavity modes for which p ^ (E 2 ), except in 
special cases. Another difference is that the response of the spin system, measured by a 
quantity analogous to our AP 2 , seems to track the energy density whenever the system is 
far from resonance, either w > O or w C fi, where uj is the frequency of the excited mode 
of the photon field and Q is the resonant frequency of the spin system. By contrast, we 
found a correlation of (E 2 ) with ^2/^2(0) only near resonance, and a correlation with AP2 
only far below resonance. This difference can be traced to the use of adiabatic switching in 
Ref. [lOj. At t = to, the spin system is assumed to be in the lower energy state, and then 
one takes the limit to ~~ * —00, and averages oscillating quantities. In this limit, for example 
sin 2 [(l/2Kfi -cu)(t- t )} -> 1/2, if Q ^ uj. This is used to go from Eq. (4.8) to Eq. (4.9) 
in Ref. [lOj. After this is done, it is not possible to let uj — > Vt. For example, the quantity 

sin 2 [(l/2)(fi - lo) (t - t )} \{t - t f ,for uj —> ft, with t fixed, asinourcase 

(D, — to) I |(f2 — u)~ 2 , adiabatic switching . 

This describes a system which was prepared in the lower state and then coupled to the 
radiation field in the distant past. In this case, one obtains expressions such as Eq. (4.12) in 
Ref. [hJ, with factors of l/(fi — uj) 2 ,1/(Q + uj) 2 , and l/(fi 2 — uj 2 ), which are proportional 
to (E 2 ) (or to the energy density, p), only if uj ^> Q or uj <C Q. 

The adiabatic switching assumption may be appropriate for a system coupled to photons 
in a plane wave mode, but does not apply to our model. The entrance of the atoms into a 
cavity is more accurately described by the sudden approximation employed here. 



IV. SOME NUMERICAL ESTIMATES 



A. Criteria to be Fulfilled 

In this section, we will analyze the feasibility of observing the suppression of the de- 
excitation probability. We consider an atom passing through a cavity with a single mode 
excited. Let us first list the desired criteria for a measurement of P 2 . We wish to impose 
the following conditions: 

(1) uj pa As. This will allow us to apply Eq. (I3"7|) , and have a situation where the de- 
excitation probability tracks the mean squared electric field. 

(2) It will be convenient to consider Rydberg atoms, which are often used in cavity QED 
experiments. Their transition frequencies are typically in the microwave range. Con- 
sider the specific case of the transition between the n — 51 and the n = 50 energy 



levels, with a transition frequency of / = 51.1 GHz and wavelength A ~ 6 mm [19 



(3) The excited cavity mode is the lowest mode, u = nJl/a 2 + 1/d 2 , so 



f U - ° / 1 1 (7ft\ 
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This requires the two longer dimensions of the cavity, a and d, each to be at least 
3 mm. 

(4) For the atom to fit in the cavity, we need the smallest dimension of the cavity, b, to 
be larger than the size of the atom. The size of a Rydberg atom with n = 50 is about 
100 nm. So we need 

b > 100 nm = 0.1 /mi. (71) 

(5) We want the travel time for the atom to traverse the cavity dimension b, to be no 
more than the time interval in which (E 2 ) < 0, which is a fraction of the period of the 
cavity mode. 

(6) To avoid the complications of relativity, and so that the atom sees a small Doppler 
shift of the cavity frequency, we want 

- < 1 . (72) 

c 

(7) We would ideally like the lifetime, r, of the excited state to be roughly equal to the 
transit time of the atom across the cavity. For the Rydberg atoms under consideration 
(n = 50), r « 3.6 x 10~ 2 sec. This would make PziO) of order unity and maximize the 
effect we are studying. 



B. Order-of-Magnitude Estimates 

Let us consider the case of the vacuum plus two photon state discussed in Sect. IHIB II 
Such a state might be created by a process which has a finite probability to emit a pair of 
photons into the lowest mode of the cavity. We have seen that in this state it is possible 
to have P^/Piiti) ~ 0.55, so there is a significant reduction in the decay rate. The duration 
of the interval of negative mean squared electric field is, from the discussion at the end of 
Sect. HUB 11 

At « — w 3 x 10- n s. (73) 
6 J 

Let us choose the smallest dimension of the cavity to be 

6^1/im, (74) 
which easily satisfies criterion (4). This requires that the speed of the atom be 

b 

v > — « 3 x 10 5 m/s. (75) 
At 1 y 1 

This can be satisfied by non-relativistic speeds compatible with criterion (6). 

The one criterion listed above which cannot be satisfied is (7), as the atomic lifetime is 
necessarily long compared to the transit time. This means that each atom has a probability 
of less than 10 -9 of decaying while transiting the cavity. However, this need not be a serious 
problem if a sufficiently large flux of atoms can be used. 
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V. SUMMARY AND DISCUSSION 



In this paper, we have considered the effects of vacuum fluctuation suppression on the de- 
excitation probabilities of atoms. This suppression occurs in quantum states of the radiation 
field in which the renormalized expectation value of the square of the electric field operator 
is periodically negative. Examples of such states include the vacuum plus two photon state 
and the single-mode squeezed vacuum state. Both of these are examples of quantum states 
exhibiting sub-vacuum effects, such as negative energy density or negative mean squared 
electric field. 

We have treated a model detector which consists of an atom in an excited state traversing 
a cavity containing photons in a non-classical state. We calculated the ratio, P 2 /P 2 (0), of 
the de-excitation probability in an arbitrary single-mode cavity field quantum state to the 
same probability in the vacuum state. This calculation assumed that the atom and field 
are near resonance, with the transition frequency approximately equal to the lowest cavity 
frequency, so that the atom had only one available decay mode. Our results showed that 
near resonance, the ratio P 2 /P 2 (0) tracks {E 2 ). Hence P 2 /P 2 (0) will have its minimum value 
when (E 2 ) is maximally negative. In certain cases, this is also when the energy density is 
maximally negative. 

We also calculated the difference, AP 2 = P 2 — p2(0), of the de-excitation probability 
between our arbitrary state and the vacuum state, where the effects of the non-excited 
modes would cancel out. In this case AP 2 = P 2 — P 2 (0) is proportional to (E 2 ) far below 
resonance, but not in other cases. In particular, only when the lowest cavity mode frequency 
was much less than the transition frequency of the atom did we find that AP 2 < when 
(E 2 ) < 0. In this case, we were able to derive a quantum inequality lower bound on AP 2 . 

Finally, we discussed the plausibility of an experiment involving Rydberg atoms to detect 
(E 2 ) < and possibly negative energy density. Although challenging, such an experiment 
may be possible. 
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APPENDIX A: BOUND ON AP 2 



Here we will derive a quantum inequality bound on the amount by which the de-excitation 
probability may be suppressed by quantum field effects. From Eq. (1361) . we have 



( 7j |P 2 (x ,t)| 7l ) = / 2 (xo) 



2(n)+2^ yj(n + 2){n + 1) Re(c n c* n+2 < 

n=0 



(Al) 



Let us examine the term in brackets. We can always choose the phases of the c's so that 
the magnitude of the second term will be largest at t — 0. Let us call the smallest possible 
value of the bracketed term S, i.e., 



S = 2(n) +2j2\/(n + 2)(n + l) Re(c n c\ +2 ) . 

n=0 



(A2) 
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Following the argument given on p. 230 of Ref.[4j, we will prove that S, and thus (E 2 ), are 
bounded from below. Since S also appears in our expressions, Eq. (|37|) and Eq. (H6l) . we can 
show that these quantities are bounded below as well. Using the fact that 



( n ) =J2 n \ C nf , (A3) 



n=0 



we may expand S as follows, 



S = 2Y, n\c n \ 2 +^^/(n + 2)(n + l) c n c* n+2 + J2\/( n + 2 )( n + 1 ) c*nC n+2 . (A4) 

n=0 n=0 n=0 



Rewriting the right-hand side of Eq. (IA3[) . we have 

oo oo 

n n l c «l 2 = H( n ~ 2 ) l c «- 



n=0 n=2 

oo 



XX n_ 1 ) K-2I - 51 l C «"2| 

ra=l n=2 

00 

-l + ^(n-l)|c n _ 2 | 2 . (A5) 



n=l 



Note also that 



£ y/(n + l)(n + 2) c n c\ +2 = ]T VT^l 41 q_ 2 c* 

n=0 1=2 

oo 



Jn(w - 1) c n _ 2 c* n . (A6) 



n=l 



To get the first equality, we let n + 2 — > to get the second equality, we relabeled £ — > n, 
and used the fact that the n — 1 term does not contribute anything to the sum. If we now 
substitute these expressions into Eq. (IA4I) . and relabel appropriately, we can write S as 



S = + [™|c„| 2 + (n- 1) |c n _ 2 | 2 + - 1) c n _ 2 c* n + yVi(n - 1) c* n _ 2 c r , 



n=l 



= -1 + 51 Iv^Cn + - lc n _ 2 | . (A7) 

n=l 

Hence 5 > — 1, and so we have that 

(E 2 ) > -/ 2 (x ) , (A8) 

and when to <C e, from Eq. (|4"Sj) . 

A ^ £ _1K^W^!W (A9) 
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